Novel metastable metallic and semiconducting germaniums 
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By means of ab initio metadynamics runs we explored the lower-pressure region of the phase 
diagram of germanium. A monoclinic germanium phase with four-membered rings, less dense than 
diamond and compressible into /3-tin phase (tI4) was found. A metallic bet- 5 phase, mechanically 
stable down to room conditions appeared between diamond and tI4. mC16 is a narrow-gap semicon- 
ductor, while bct-5 is metallic and potentially still superconducting in the very low pressure range. 
This finding may help resolving outstanding experimental issues. 



INTRODUCTION 

The fundamental character and technological relevance 
of group-IVa elements (tetrels) have motivated repeated 
investigations and systematics on their polymorphism p]- 
[5] . Carbon polymorphs are promising as hard and trans- 
parent materials [6j [7]. Semiconducting silicon (Si) is 
versatile for micro- and nanoelectronic devices, while ger- 
manium (Ge) displays comparatively higher carrier mo- 
bility, a finer band gap tunability, and good compatibility 
with high-dielectric constant materials |8] . Metallization 
occurs in silicon and germanium upon compression pQ. 
In Ge lowering of phonon frequencies promotes electron- 
phonon coupling towards superconductivity [2j [3] . The 
possibility of metallic germanium under room conditions 
is very intriguing and intensively debated [2j [9], while 
superconductivity in elemental Ge appears under pres- 
sure [3 J. In the lower pressure range, improved optical 
properties via band-gap tuning can be achieved in a dif- 
ferent polymorph. 

Engineering viable new compounds with superior prop- 
erties entails a detailed understanding of structural 
changes [10 j. Under pressure germanium bears simi- 
larities with silicon pQ by comparatively higher transi- 
tion pressures with respect to Si, due to Ge core d- 
electrons [11 . Upon compression semiconducting Ge (cu- 
bic diamond) transforms into /3-t'm type (space group 
I4i/amd) at about 10 GPa [12] , and then to Imma 
phase [13] , simple hexagonal (P6/mmm) [14 , followed 
by orthorhombic Cmca phase [15] and finally upon fur- 
ther compression above 180 GPa, by the hexagonal close- 
packed arrangement (P6s/mmc) [15] . 

The phase diagram of germanium is further compli- 
cated by a family of tetrahedral structures p2][T7]- Type- 
II clathrate Ge(cF136) exists at ambient conditions [17] . 
Other germanium modifications are reported, Ge(tP12), 
Ge(cI16) (7-silicon type, BC8) and Ge(hR8) [16J. BC8 
Ge [18] is accessible through decompression from /3-tin 
Ge. In a nutshell, during the last few years, new dense 
and open phases of germanium have been experimentally 
observed or synthesized. Nevertheless, a systematic ap- 
proach to including known and finding novel germanium 




FIG. 1. Lower pressure region of the phase diagram of Ge, 
augmented by two novel phases mC16 and bct-5, found by ab 
initio metadynamics runs, bct-5 shows characterisitic square 
pyramidal 5-fold coordination of Ge atoms. In monoclinic 
mC16 four-rings are a characteristic feature. The arrows indi- 
cate the direction of metadynamics evolution. The pressures 
were evaluated based on the common tangent construction 
(see below, Fig. 



forms is still outstanding. 

In this Letter we explore the energy landscape of ger- 
manium, both at ambient conditions and upon moder- 
ate compression. By means of ab initio metadynamics, 
we efficiently sampled structural transformations along 
appropriate collective reaction coordinates. Besides the 
already known dense phases of germanium, we found two 
novel allotropes (green in Fig.[l]). The first one is a mon- 
oclinic modification of germanium (mC16) slightly less 
dense than diamond, is an indirect band gap semiconduc- 
tor and is unprecedented for tetrel elements. The second 
one is a five-coordinated (square pyramidal) metallic in- 
termediate structure (tI4, sp. gr. I4/mmm), which in- 
curs in the diamond (cF8) —> /3-tin phase (tI4) transition, 
and which has been postulated to exist in homologue sil- 
icon (bct-5) [19 . We report on transformation paths, 
energetic, mechanical and electronic properties. For the 
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metallic bct-5 phase, we calculated superconducting tem- 
peratures down to ambient pressure compared to /3-tin 
phase, based on the electron-phonon coupling mecha- 
nism [20] . 

METHODS 

Different approaches are successful in the prediction 
of novel polymorphs of the elements [2TJ [22] . Important 
discoveries have been achieved by means of random tech- 
niques, genetic (evolutionary) algorithms, or accelerated 
molecular dynamics [10] . Metadynamics [23-25 allows 
for the exploration of the energy surface along one or 
more collective reaction coordinates. The method is in- 
dependent of the level of theory used, it does not require 
prior knowledge of the energy landscape and its sam- 
pling efficiency can be enhanced by parallel runs started 
from different configurations. The time-evolution of the 
system is biased by a history-dependent potential con- 
structed as a sum of Gaussians centered along the tra- 
jectory, in order to discourage the system from visiting 
already harvested regions of the potential [26 . 

All metadynamics runs were performed with at least 
eight atoms in the simulation box which served as a col- 
lective (6-dimensional) variable. This minimal box ap- 
proach was successfully used in the prediction of novel 
carbon polymorphs, recently [7 . The size of the mini- 
mal box ensured commensurability of all already known 
phases (except for the clathrate II phase, which requires a 
minimum of 34 atoms) either open or dense. Each meta- 
dynamics metastep consisted of molecular dynamics runs 
in the NVT ensemble for 0.5 ps (timestep 2 fs) at 300 K. 

Metadynamics was performed with different molecu- 
lar dynamics layers. A Density Functional Tight Bind- 
ing (DFTB) [27] level of theory, as implemented in the 
T-point-only DFTB module of the CP2K code [28l [29], 
ensured rapid and accurate sampling in the low-pressure 
regime, characterized by four-connected Ge atoms. For 
higher-pressures SIESTA [30 was used as the DFT 
molecular dynamics layer, allowing for k-point runs. 
Electronic states were expanded by a single-^ basis set 
constituted of numerical orbitals with a norm-conserving 
Troullier-Martins [3Tj pseudopotential description of the 
core levels. Single-^ basis set dramatically reduces com- 
putational times providing nonetheless, the right topol- 
ogy and energy differences of all the Ge allotropes under 
study. The charge density was represented on a real- 
space grid with an energy cutoff [30] of 200 Ry. A 
Monkhorst-Pack k-point mesh of 2 x 2 x 2 ensured the 
convergence of the electronic part. High-pressure meta- 
dynamics was performed based on DFT. Lower pressure 
regions were initially explored by DFTB, followed by 
DFT metadynamics upon discovery of interesting novel 
polymorphs. In the lower pressure range the transferabil- 
ity between DFTB and DFT is unflawed. 



Electronic structure, phonon dispersion curves and su- 
perconducting properties were calculated with the Quan- 
tum Espresso (QE) [20] [32] package. The superconduct- 
ing critical temperature T c was evaluated based on the 
Allen and Dynes modification of the McMillan formula. 
This required calculating the electron-phonon coupling 
strength A via the Eliashberg function. The Coulomb 
potential value was /i=0.1. A q-mesh of 8 x 8 x 8 was 
used for the evaluation of the dynamical matrix, while a 
Monkhorst-Pack k-point mesh of 12 x 12 x 12 ensured 
convergence of the electronic part. 

The structures visited during each run were character- 
ized by their vertex symbols, which contain the informa- 
tion on all the shortest rings meeting at each atom, and 
coordination sequences, as implemented in the TOPOS 
package [33] . In case of new structures ideal space group 
and asymmetric units were identified with the Gavrog 
Systre package [34]. Subsequently a variable-cell ge- 
ometry optimization was performed (DFT-GGA, PBE 
functional [35]) in a plane- wave pseudopotential frame- 
work [32] using Vanderbilt ultrasoft pseudopotential as 
supplied by Perdew-Zunger with non-linear core correc- 
tion [36] [37] . A k-point mesh of 8 x 8 x 8 ensured con- 
vergence of the electronic part, while a plane- wave basis 
set with an energy cut-off of 30 Ry was applied. 

RESULTS AND DISCUSSION 

The mC16 structure (Fig. [TJ C2/m, a=7.6094 A, 
6=7.9746 A, c=6.5668 A, /3=104.10°) arised from a meta- 
dynamics run started from diamond (8 atoms box, p=l 
bar, T— 300 K). Ge atoms occupy three Wyckoff posi- 
tions: (4z) 0.70984 0.50000 0.67434, (4z) 0.60981 0.50000 
0.29080, (8j) 0.65012 0.76388 0.11596. Strikingly, mC16 
is less dense than diamond (see Fig. [2|, although topo- 
logically as dense as diamond or lonsdaleite. Its bulk 
modulus amounts to 51.2 GPa that is slightly lower than 
that of the diamond type- structure (60.7 GPa), estimated 
from the fit to the third order Birch-Murnaghan equa- 
tion of state. Applying pressure to the mC16 allotrope 
in an additional metadynamics run resulted into a direct 
transition to the /3-tin phase. Decompressing the latter 
is known to generate metastable phases typically denser 
than diamond Ge. Therefore, a viable route to mC16, 
like for other recent germanium allotropes, could rather 
be the oxidation of suitable germanium Zintl salt precur- 
sors, i.e. via chemical synthesis. 

Upon compression diamond transforms into /3-tin and 
it subsequently follows the same transition sequence of 
silicon phases. Along the diamond —> /3-tin transition 
metadynamics (64 atoms box, p= 10 GPa, T— 300 K) 
visited an intermediate of bct-5 (sqp) topology (I4/mmm, 
a=3.5491, c=6.4478, Ge(4e) 0.0 0.0 0.19273). The bct- 
5 bulk modulus is 58.7 GPa, slightly lower than that of 
the /3-tin phase (68.2 GPa). This five-connected struc- 
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FIG. 2. Equation of states of Ge mC16 and bct-5, compared 
to cubic diamond and /3-tin type, bct-5 features a reduced 
volume per atom compared to diamond type (cF8) , while the 
total energy minimum lies lower than /3-tin. mC16 on the 
contrary is less dense and energetically close to the diamond 
type. The tangents used for evaluating equilibrium pressures 
of Fig. [I] are highlighted. 

ture has been proposed for silicon [T9] . but has never 
been observed so far. The total energy /volume curves of 
Fig. [2] suggest bct5 as a conventional product of diamond 
compression. However, under hydrostatic conditions /3- 
tin is formed from diamond, while decompression leads 
to other germaniums, although indications of minority 
phases exist. 

The transformation affects only one box parameter, 
suggesting nonhydrostatic shearing as the protocol of 
choice towards bct-5, also supported by the magic-stress 
approach which led to bct-5 in silicon. Alternatively, low- 
temperature compression may be considered. The evolu- 
tion of the enthalpy profile of metadynamics runs from 
diamond Ge (Fig. [4| shows bct-5 as a narrow plateau 
around metastep 160. Metadynamics runs started from 
bct-5 (Fig. pi black curve) confirms it as a proper in- 




FIG. 3. Evolution of the bonding situation from diamond to 
/3-tin type. The ELF map is showing four bond attractors 
for diamond Ge (a, 77 = 0.58), one + four bond attractors for 
bct-5 (b, 77 = 0.53, transparent green isosurface 77 = 0.48), 
two + four bond attractors for /3-tin (c, 77 = 0.51). 
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FIG. 4. Metadynamics (DFT) runs in a 64 Ge atoms box 
(pi = 10.0 GPa, T=300 K). bct-5 appears as a stepwise feature 
in the enthalpy profile of the run started from diamon (red 
line). Runs commenced from bct-5 (black line) evolve into 
/3-tin. Configurations corresponding to distinct points along 
the runs are detailed below the graph. 



termediate along the transition, which can be quenched 
down to room pressure, and which is mechanically stable 
(see below). Mechanistically, the coordination number 
increases from 4 to 5 on shortening one bond, followed 
by flattening of the pristine tetrahedron and formation 
of the square pyramidal (sqp) geometry of bct-5, Fig. [4)3. 
The four bonds in the pyramid basis are 2.62 A long, the 
axial 2.48 A. 

The evolution of the bonding situation from Ge dia- 
mond to /3-tin over bct-5 is shown in Fig. [3] Calculation 
of the ELF [38] shows four, one + four and two + four 
bond attractors, respectively. The five "bonds" in this 
orbital-deficient, electron-deficient metallic bct-5 result 
from the sp Ge valence shell. This bonding scenario is 
reminiscent of the recently discovered superconducting 
Zintl phase CaGe3 [39], isosymmetric with bct-5. 

The electronic band structures and phonon dispersions 
of mC16 and bct-5 are shown in Fig. [5] The tetrahedral 
phase is semiconducting with an indirect band gap of 1.43 
eV (PBE-GGA [35]), while bct-5 is metallic and stable 
down to GPa. Isothermic-Isobaric molecular dynamics 
(1 bar, 300 K, 2.5 ps) confirmed the stability of bct-5. 
mC16 is characteristic due to the presence of four-rings. 
However, this does not imply overall structure destabi- 
lization [40]. The expectation of a strained geometry is 
in fact not supported by total energy calculations, which 
place mC16 among the energetically lowest Ge allotropes. 
The indirect band gap and the low density (compared to 
the diamond type) makes this germanium an attractive 
material. The need for a "negative" pressure makes a 
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FIG. 5. Band structures and phonon spectra [32] (0.0 GPa) 
of mC16 (left) and bct-5. mC16 is a narrow-gap semiconduc- 
tor (band gap = 1.43 eV), while bct-5 is metallic. Both are 
mechanically stable. 



FIG. 6. Evolution of u, A and T c as a function of pressure 
for bct-5 and /3-tin, calculated based on the electron-phonon 
coupling model. The calculated equilibrium pressure between 
bct-5 and /3-tin marks the boundary between the phases. The 
model predicts an increase of T c in the lower pressure region. 



chemical path plausible. 

The prominent property of bct-5 is the conservation 
of metallic character down to ambient conditions. Cal- 
culations and experiments have shown an increase of the 
superconducting temperature on lowering pressure, with 
superconductivity still present around 8,7 Gpa (T c = 5 
K) [2]. The evolution of T c as a function of pressure 
for bct-5 is shown in Fig. [6] The ongoing debate on 
the possibility of "strange metallic" states other than 
/5-tin [2j |9], and the need for further explanations of 
the survival upon decompression of certain "metallic" 
metastable states that could not be reliably assigned to 
any known Ge phase, open the door for a serious consid- 
eration of the role of the bct-5 in the lower pressure range, 
as a metallic state that qualifies for higher T c values (via 
the McMillan relation). 

In conclusion, with metadynamics runs we have found 
and characterized two novel Ge polymorphs. The first 
one, mC16, is an indirect gap semiconductor, and its 
structure is unprecedented for the tetrel family. The sec- 
ond bct-5 polymorph was suggested for Si. Our simula- 
tions lean strong relevance to bct-5 in the lower pressure 
range, as a further metallic superconducting phase ca- 
pable of stability at room conditions. We expect our 
predictions to stimulate further experimental work. 
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